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A frustrated Ising model on a diamond hierarchical lattice is studied. We obtain the exact 
partition function of this model and calculate the transition temperature, specific heat, entropy, 
magnetization, and ferromagnetic correlation function. Depending on the magnitude of a pa- 
rameter giving the frustration, there exist three types of ground states: ferromagnetic, classical 
spin-liquid with highly developed short-range order, and paramagnetic. The dependence of the 
zero-temperature entropy on the frustration parameter has an infinite number of steps. The 
temperature dependence of the specific heat exhibits many peaks with decreasing temperature 
and entropy loss. The dominant spin configurations at low temperatures are also specified. 
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1. Introduction 

Many magnetic systems with competing interactions 
exhibit frustration that leads to multiple ground states 
called spin glass or spin liquid even at the classical level. 
These interesting properties are attributable to the del- 
icate balance of frustrated spin-spin interactions. If we 
can obtain the exact partition function of a model with 
frustrated magnetic phenomena, we can obtain a better 
understanding of complicated frustrated phenomena. 

In statistical physics, few exactly soluble models with 
phase transitions are known, such as two-dimensional 
Ising model, 1 eight-vertex model, 2 one-dimensional van 
der Waals gas model, 3 and a hierarchical model. 4 Since 
Berker and Ostlund proposed a hierarchical model re- 
lated to the renormalization group method, ■ 5 many dif- 
ferent models on hierarchical lattices have been proposed 
and developed. 5-8 The magnetic and thermodynamic be- 
haviors have been studied, 9-11 and the distribution of 
the zeros of the partition function and the critical ex- 
ponents have also been obtained. 10 ' 12 Using a hierarchi- 
cal model with competing ferro- and antiferromagnetic 
interactions, McKay et al. studied the spin-glass behav- 
ior 13 and Nogueria et al. investigated the local magneti- 
zation. 14 However, to the best of our knowledge, no hier- 
archical model describes the spin-liquid ground state. In 
this study, we consider a hierarchical Ising model with 
frustrated interactions that lead to the classical spin- 
liquid ground state. 

This paper is organized as follows. We introduce our 
diamond hierarchical lattice and frustrated Ising model 
in §2 and describe our recursion relations in §3. Ther- 
modynamic and ground state properties are described in 
§4 and §5, respectively. In §6, the temperature evolution 
of dominant spin configurations at low temperatures is 
discussed. In §7, we summarize the results obtained in 
this study. 

2. Lattice and Hamiltonian 

Our diamond hierarchical lattice is constructed by the 
infinite iteration procedure shown in Fig. 1. We begin 
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1st stage 2nd stage 



(n-l)th stage 





Fig. 1. Construction of the hierarchical lattice. The n-stage lat- 
tice comprises four (n — l)-stagc lattices represented by the gray 
diamonds. The sites denoted by A, B, 1, and 2 are contact sites 
between two of the (n — l)-stage lattices. 



with a single bond expressed by the solid line, and call 
this the first stage. We replace this single bond by a 
diamond-shaped basic unit comprising four solid lines 
and one dotted line to obtain the second-stage lattice. In 
order to obtain the third-stage lattice, each of the solid 
lines in the second-stage lattice is replaced by a diamond 
unit, and the dotted line is left untouched. In general, 
the n-stage lattice is constructed by replacing each of 
the solid lines in the (n — l)-stage lattice a the diamond 
unit. The dotted lines are always left untouched when 
we proceed to the next stage. It should be noted that 
the n-stage lattice comprises four (n — l)-stage lattices. 
The n-stage lattice has four contact sites between the 
(n — l)-stage lattices, and these are indicated by A, B, 
1, and 2 in Fig. 1. We call sites A and B (1 and 2) as 
surface (internal) contact sites. The n-stage lattice has 
N n = |(4™- 1 + 2) sites, 4"" 1 solid lines, and i(4"- 1 - 1) 
dotted lines. 

In order to define a frustrated Ising model on the di- 
amond hierarchical lattice, we place Ising spins on each 
of the lattice sites. We call spins a a and gb defined on 
the surface contact sites as "surface spins," and spins cti 
and (72 defined on the internal contact sites as "inter- 
nal contact spins." Spins on both ends of a solid and a 
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dotted line couple fcrromagnctically and antifcrromag- 
netically by — J and aJ, respectively. The Hamiltonian 
can be written as 

H = - J (JiOj + aJ UiGj - H crj, (1) 

where the first and second sums run over all pairs on the 
solid and dotted bonds, respectively. Although we study 
the zero- held properties, we have introduced the Zccman 
term in our Hamiltonian to calculate the magnetization 
and correlation functions. 

Antifcrromagnetic coupling tends to destabilize the 
ferromagnetism induced by the ferromagnetic coupling. 
At finite temperatures, there exists a continuous phase 
transition between the ferromagnetic and disordered 
phases as a function of a. At zero temperature, the phase 
transition cannot be continuous because no fluctuations 
exist. In this study, we show that an infinite number of 
first-order phase transitions occur at zero temperature 
instead of a continuous phase transition at finite temper- 
atures. Disordered ground states near the ferromagnetic 
transition point are regarded as spin-liquid because the 
ferromagnetic short-range order is highly developed. The 
residual entropy is step-like as a function of a, leading 
to unique thermodynamic properties. 

3. Recursion Relations 

Here, we derive recursion relations for our hierarchical 
model. The partition function of the n-stage lattice is 
given as 



Paper 
where 
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Z n = Qn(<TA,<J B )e h{aA+aB \ 



(2) 



where 



Qn{?A,<T B ) = A n e K ^ A ° B+h ^ aA+a ^ (3) 

and h = H/T. (We use the unit with k& = 1 in this 
paper.) For n = 1, we have Ay = 1, Ky = K (= J/T), 
and hi = 0. For n > 2, noting that the n-stage lat- 
tice comprises four (n — l)-stage lattices, we can express 
Qn{&A, fs) in terms of the (n— l)-stage lattice as follows: 

Qn(<7A, 0\b) = Qn-l{&A, &l)Qn-l{<Tl, & b) 

(Tl,(T2 

xQn-l(&A, cr 2 )Q„_i (cr 2 , Ub) 

v e -aKcr 1 a2+h(<j 1 +a 2 ) 

= 2A*_ 1 {e aK + e~ aK 

x cosh[2A„_i {a A + a B ) + ih n -y + 2h}}. (4) 

Using eqs. (3) and (4) , we obtain the following recursion 
relations: 

K n = f K {K n -y,h n -y,h) + fK(K n -l, -hn-i,—h) 

-2f K {0,hn-uh), (5) 

h n = 2h n -i + f K (K n _i,h n -i, h) 

-fK(K n -i,-h n -i,-h), (6) 

and 



f K (K,h,h) = - ln[e aK + e.- aK cosh(4K + 4h+2h)]. (8) 

We express Ki , hi, and In Ai = Ai as a power series in 
h as follows: 

K^Kf +K? ) h 2 + 0{h% 



hi = h^'h + 0{h A ), 

I^lf + Af ) h 2 + o(h 4 ). 

The recursion relation for is given by 



cosh(4ifW 1 ) 



K<® =-ln- 

2 2cosh(aA) 



(9) 
(10) 

(11) 
(12) 



We can write A^ and h n ^ in terms of as follows: 



(0) 



4 n-l _ 1 



ln[4cosh(aA:)], (13) 



n— 1 n—1 



4 1)= $>* ; ; 2(i+ Ps )+ Pn _ 1 , (u) 



i— 1 s—i+1 



where 



e~ aK sinh(4fiT. (0) ) 



Pi 



(15) 



e aK + e- aK cosh(4K<> 0) )' 
In a similar manner, we have 

n—2 n—1 

*£°=I> II 2 Ps + Qn-l (16) 



i— 1 s—i+1 



with 



J tanh(aA') - 1 e 2aK cosh(4A; (0) ) + 1 1 
* ~\ 2 [cosh(4iff ) ) + e 2aK } 2 J 

x(l + 2ftf) 2 , (17) 



and 



A {2) = ^ 4 »-{^-(l+/ l «) 2 [tanh(a^)-l]}. (18) 

i=2 

We use K„, K n 2 \ h n , A„, and A^ obtained from 
eqs. (12)-(18) to study the phase diagram, specific heat, 
entropy, magnetization, and ferromagnetic correlation 
function in the next section. When calculating these 
quantities, we choose the maximum value of n as 20~25 
and check if the calculated results can be regarded as 
those in the thermodynamic limit. 



\nA n =41n^„_i+A'„ + 2/^(0, ft), (7) 
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4. Thermodynamic Properties 

4-1 Phase diagram 

The transition temperature and phase diagram can be 
determined from the fixed point of cq. (f2). The phase 
boundary in the a — T plane is given by 



(a) 




f 



(19) 



The phase diagram is shown in Fig. 2. When a < 1, a 
second-order transition from the ferromagnetic phase to 
the disordered phase occurs at finite temperature. For 
1 < a < 2, it will be shown later that the ground state 
becomes a classical spin-liquid state in which the short- 
range order is highly developed but the magnetization is 
zero and the entropy is finite. When a > 2, the ground 
state is paramagnetic in which the short-range order van- 
ishes. The paramagnetic ground state consists of inde- 
pendent dimcrs composed by the coupling aJ. 



disordered phase 




Fig. 2. Phase diagram in the a — T plane. At zero temperature, 
a short-range ferromagnetic order develops for 1 < a < 2 (spin- 
liquid ground state) but not for 2 < a (paramagnetic ground 
state). 
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4-2 Specific heat and entropy 
Using the relation 

lim Z n = 4A£ 0) cosh K 
we obtain the specific heat per spin 



(20) 



C= lim 



K 2 d 2 



N n dK 2 
and the entropy per spin 



-[il )+ln(cosh^ ))], (21) 



,5' 



lim — — 

n— >oc JV„ 



l-^)[il 0) +ln(co S hif(°))]. (22) 



dK 



Using the recursion relations for Kn arid An' m 
cqs. (12) and (13), we can evaluate C and S at finite 
temperatures. 

The dependence of the specific heat on temperature 
is shown in Fig. 3. The results for a = 0~1, for which 
a ferromagnetic phase appears at low temperatures, are 
shown in Fig. 3(a). The specific heat for a = diverges at 
the ferromagnetic transition temperature T c = 1.641J. 9 



AO) 



Fig. 3. Temperature dependence of the specific heat for (a) < 
a < 1, (b) 1 < a < 2, and (c) a > 2. In (b), the original 
curves are shifted upwards by multiples of 0.1, and the number 
corresponding to each curve represents the value of a. 



When a increases, the divergence point decreases, which 
agrees with the behavior of the above-obtained phase 
boundary. For a = 0.8, we find a broad maximum at 
T = 0.66J, and the ferromagnetic transition tempera- 
ture T c = 0.5545J is lower than the broad maximum 
temperature. For a = 0.9, many peaks appear at tem- 
peratures slightly greater than the transition tempera- 
ture T c = 0.2881J. Figure 3(b) shows the specific heat 
for a = 1 ~ 2. We find that the result for a = 1 has 
an infinite number of broad peaks. For a > 1, the num- 
ber of broad peaks becomes finite. The evolution of the 
specific heat curve as a function of a exhibits a complex 
behavior, and it is described in detail in relation to the 
entropy later. Figure 3(c) shows the results for a > 2. In 
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Fig. 4. Contour map of the specific heat in the a — T plane for 
a > f. 



this figure, we choose a J as the unit of temperature. The 
result for a = 2.3 has two peaks. When a is increased, the 
low-temperature peak merges with the high-temperature 
one. Then, with a further increase in a, the curve moves 
toward the a = oo limit monotonically. In Fig. 4, we 
show the results shown in Figs. 3 (b) and (c) in the form 
of a contour map in the a — T plane. 

It is helpful to calculate the entropy in order to under- 
stand the above mentioned behavior of the specific heat. 
The entropy per spin as a function of a at low temper- 
atures is shown in Fig. 5. We have S = in the ferro- 
magnetic ground state for small a and S = (1/2) log 2 ~ 
0.347 in the paramegnetic ground state for large a at 
T — > 0. We find that many steps appear between them. 
We define the positions of steps at T — ► as ot\, a%, ■ ■ ■ 
in descending order. (In the next section, it is shown that 
ai = 2 l /(2 l — 1).) When a decreases and passes a/, an 
additional peak appears in the temperature dependence 
of the specific heat, as seen in Fig. 3(b). Figure 4 clearly 
indicates that we can ascribe each low-temperature peak 
in the specific heat to a step in the a-dependence of en- 
tropy. For example, in Fig. 3(b), the peaks (or shoulders) 
of a = 1.9, 1.7, 1.5 indicated by the open circles corre- 
spond to the entropy step at a = oti (= 2), and those of 
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a = 1.3, 1.2 indicated by the closed circles correspond to 
the entropy step at a — ct2 (= 4/3). The entropy exhibits 
a spike at around a = 2. as seen in Fig. 5. Then, the spe- 
cific heat for a slightly greater value of a may have an ad- 
ditional peak. In fact, the specific heat curve for a = 2.3 
shown in Fig. 3(c) has an additional low-temperature 
peak. The next and more important problem is why the 
a-dependence of entropy has such a structure. This issue 
is discussed in the following section. 

4-3 Magnetization 

Griffiths and Kaufman proved that the hierarchical 
model has a thermodynamic limit. 6 The surface spin 
state does not have a significant effect on the physical 
properties of the system in the limit of n — ► 00. There- 
fore, instead of an extremely small external field applied 
on this system, we use a a = 1 and as = 1. The magne- 
tization per site, m, is defined by following expression: 
1 d 

lnQ n (a A = l,a B = 1) 



lim 



h->0 ah 



= lim 



3h 



(i) 



If the system does not have a long-range order, no magne- 
tization appears under this setting of boundary surface, 
and if it does, a finite magnetization can be obtained in 
the limit of a vanishing external magnetic field. 

The temperature dependence of magnetization for < 
a < 1 is shown in Fig. 6, where the horizontal axis 
represents the reduced temperature, (T c — T)/T c . For 
a = 0, we obtain a magnetization index [3 = 0.1617, 
which agrees with the result obtained by Margado et al. 9 
With a decrease in a, the form of the magnetization curve 
tends to resemble a step function, and the index /3 de- 
creases and vanishes for a — > 1 (see the inset in Fig. 6). 
The continuous variation of the critical exponents has 
been reported in previous studies. 6 ' 7 




0.0 0.2 0.4 0.6 0.8 1.0 



a = 0.0 
a = 0.4 
a = 0.7 
a = 0.8 
a = 0.9 



Fig. 5. Entropy 5 as a function of a at T = 0.1J and 0.01J 



0.6 0.4 

(T c -T)/T c 



Fig. 6. Magnetization as a functions of temperature. The inset 
shows the o-dependence of the critical index /3, defined by m k 
(T c — T) 13 around T = T C . 



4-4 Ferromagnetic correlation function 

We calculate a ferromagnetic correlation function to 
study how the short-range order develops in the disor- 



J. Phys. Soc. Jpn. 



Full Paper 



H. Kobayashi, Y. Fukumoto and A. Oguchi 5 



dcrcd region. Using 



(24) 



we define the ferromagnetic correlation function as 
(A|) = lim (A|(n)) 

r 1 ,. d 2 \nZ n 
= lim — — lim 



n->oo N n /i-»o dh 2 
3 



= lim 



i(2) + K (2) 



1 



(25) 



In Fig. 7, we show the a-dependence of the correlation 
function at T = 0.01 J, 0.1 J, and J. A system with a > 1 
is in the disordered phase at all temperatures. However, 
if a ~ 1, the short-range ferromagnetic correlation is well 
developed at low temperatures. We have limT^o(Af) > 
for 1 < a < 2 and lim r _, (A|) = for a > 2. 

A distinctive feature of the present model is that a 
short-range order develops for 1 < a < 2 between the fer- 
romagnetic and paramagnetic ground states. If the fer- 
romagnetic state for a < 1 is comparable to the solid 
state of matter and the paramagnetic state for a > 2 to 
the gas state, this state for 1 < a < 2 is comparable to 
the liquid state. Therefore, we can consider this ground 
state in the region of 1 < a < 2 as a classical spin-liquid 
ground state. 



< 

V 



I 

7 = 0.01.7 
T = 0AJ 
T = J 



5.1 Ground state energy and spin configurations 

We begin with the case n = 2, and denote the 
spin configuration as (cta, <?b, o~i, a 2 )- When a is small, 
the lowest-energy spin configurations are ferromagnetic 
states, (t, Tj Ti T) and (4,1, 1, 4), and the ground state en- 
ergy is given by —4 + a. When a increases, the ferromag- 
netic states are destabilized against eightfold-degenerate 
spin configurations, {(cr^, a B , <Ji, o- 2 )\oa = X,\\°b =T,i 
; °"i =T;I; CJ 2 = ci}, whose energy is —a. The energy 
crossing point is a = 2 (= ai). Therefore, the ground 
state energy of the second-stage lattice, ^(a), is given 
by 



E 2 {a) 



—a for ot\ < a 

— 4 + a for otherwise 



(26) 



Thus, the ground state for n = 2 is divided into two sec- 
tors depending on a. We call ol\ < a as the first sector 
and the other one as the ferromagnetic sector. In addi- 
tion, the ground state energy of the first (ferromagnetic) 
sector of the second-stage lattice is denoted by E^i = —a 
{E 2J = -4 + a). 

For the third-stage lattice, there are two crossing 
points in the ground state energy E^ipi) as a function 
of a (see Fig. 8), and thus, we have three sectors: the 
first sector is a± < a; the second, a.?, < a < cti, where 
oi2 = 4/3; and the third, a < a 2 , is the ferromagnetic 
sector . The energy of the first sector, E^i = —5a, satis- 
fies £3,1 = 4i?2,i — a, which shows that each of the four 
diamonds in the third-stage lattice gives energy £2,1 and 
the internal contact spins o\ and a 2 couple antiferro- 
magnctically. For the energy of the ferromagnetic sector, 
E 3 j = — 16 + 5a, we have E 3 j — iE 2 j + a. Note that 
the internal contact spins in the ferromagnetic sector lead 
to energy loss by a. This energy loss is the cause of the 
appearance of the second sector in the third-stage lat- 
tice. For the energy of the second sector, £3.2 = —8 — a, 
the relationship £ 3 ^ 2 = 2£ 2 ,i + 2£ 2 j — a holds. This 
relationship indicates that two of the four diamonds give 
up being in the ferromagnetic state to obtain an energy 
gain from the antiferromagnetic interaction between o~\ 
and a 2 . 

For the fourth-stage lattice, we have four sectors, as 
seen in Fig. 8. It is easy to check that recursion relations 



Fig. 7. Ferromagnetic correlation function as a function of a at 
T = J, 0.1J, and 0.01J. 



5. Ground State Properties 

In the previous section, we have shown that a spin- 
liquid ground state appears between the ferromagnetic 
and paramagnetic ground states on the basis of the cal- 
culated results of the entropy and the correlation func- 
tion. Here, we study ground-state spin configurations in 
order to understand how the short-range order develops 
in the spin- liquid ground states. 




1.0 / 1.2 / 1.4 1.6 



Fig. 8. Dependence of the ground state energy per site on a for 
n = 2, 3, 4. 
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for energies similar to the third-stage lattice also hold. In 
general, we have n sectors in the ground state of the n- 
stage lattice. We denote the ground state energy of each 
sector by E n ^\ for ai < a, E„ t 2 for a 2 < a < a%, • • • , 
E n , n -i for a„_i < a < a n _ 2 , and E„j for a < a n -i- 
It is expected that the (n — l)-th sector comprises the 
(n — 2)-th and ferromagnetic sectors in the (n — l)-stage 
lattice, and each of the other sectors comprises sectors 
from the same number in the (n — l)-stage lattice; this 
is confirmed in the following subsections by comparing 
the zero-temperature entropy and correlation function 
with those obtained by the finite-temperature calculation 
presented in the previous section. 

Here, we derive expressions of the ground state energy 
for the 77,-stage lattice on the basis of the above observa- 
tion. For the ferromagnetic sector, we have 



(27) 



This recursion relation, along with the initial condition 



E- 



2 J 



-4 + a, gives 



E nJ = -4"- 1 + -(4"- 1 - 1). (28) 
For sector I = n — 1, using the recursion relation 

E n , n -i — 2E n -ij + 2i? n _i,„_2 — ol (29) 
and the initial condition E 2 ,\ = —a, we obtain 



rp on 

^nji — 1 ^ 



-(4 



n-l 



3 -2" + 5). (30) 



Finally, for sectors I = 1, • • • , n— 2, the recursion relation 
E n j = 4£ n _i,, - a (31) 



gives 



E,, 



->2n-Z-l 



3 



+ 4 n -' + l), (32) 



where we have used the expression of Ei+u obtained 
from eq. (30). The condition £",u+i = E n j yields the 
lower bound ai in the form 

2 l 

ai = 2 1- T (33) 

for 1= 1,2,- •• ,n- 1. 

In order to study the lowest-energy spin configurations 
in greater detail, it is convenient to introduce the follow- 
ing notations. We represent a spin configuration of the 
n-stage lattice as {a a, o~b,s), where a a and ctb are sur- 
face spins and s represents the configuration of the other 
spins. We define C^'" 3 ' as a set of internal configu- 
rations s that gives the lowest energy for the Ith sector 
under the condition that the surface spin state is fixed 
as (<ja,&b)- For example, setting s = (eri,cT2) for the 
second stage (n = 2), we have 



ma), am 



for the first sector (I = 1) and 



°2,/ 



{(CTA,cr B )} for OA = 

for a a 7^ o\B 



(34) 



(35) 



for the ferromagnetic sector. A detailed discussion of 



C 



(<7a,0-b) 
,1 



can be found in Appendix A. In particular, the 



derivation of an analytical expression of |C„ ; 4 <TS ^I: the 

total number of elements in the set C^^ ,aB \ is presented 
there. 



5. 2 Entropy 

The zero-temperature entropy per site for the Ith sec- 
tor at n — > 00 can be written as 

\<TA,&B / 



3- 2 1 ' 1 - 1 



log 2, 



(36) 



where we have used the expression for \C^ T f" aB \ given 
in eq. (A- 11) of Appendix A. Denoting the entropy as a 
function of a by S'(a), we have S(a) = Sooj for ai < a < 
ai-i- With regard to the entropy at a boundary a = ai, 
noting that <S'(q;) at T = 0.01J exhibits at around a = ct\ 
(Fig. 5), we expect that S{a{) > at T = 0. As 

derived in Appendix B, S(ai) can be written as 



S(ai) = SryoJ + 



8 -A 1 



log 



1 + 4(2~ 2 +4~ 2 ) . (37) 



In Fig. 9, we show the entropy obtained by the above 
expressions, along with that obtained by the finite- 
temperature calculation presented in the previous sec- 
tion. We find that both results are consistent with each 
other. A clear spike is obtained at around a = Qfi(= 2). 
The height of this spike is S(ai) — Soo.i = ^ log | — 

0. 076. There is a spike at around a = ai for I > 2; the 
height of this spike decreases rapidly with an increase in 

1. The inset shows a double logarithmic plot that shows 
the behavior in the immediate vicinity of a = 1. We find 
that the discrepancy between the results of T j J = 10~ 2 
and T/J = is pronounced below a — 1 ~ 10~ 2 . How- 
ever, the result of T/J = 10~ 3 indicates that lowering of 
the temperature resolves the discrepancy, as expected. 

5.3 Magnetization 

In order to obtain the magnetization of the ground 
state we consider the surface spins to be upward, as be- 
fore. We denote the magnetization of spin configurations 
in C ( Jf ] as M(C^' T) ). Then, we find M{C { 2 ]{ ]] ) = and 

MfC^f) = 2, and obtain the following recursion rela- 
tions: 



M{cli] ) ) = 2 + AM{C\11 J ) 



and 



M(C^ ) )=2M(C«:\\ l ) 



where C^l n _ x = C^ f . We have used eqs. (A-l), 
(A- 2), and (A- 7) of Appendix A in the derivation of these 
relations. Solving these equations, we have 



for I < n, 



(38) 



(39) 



■(4 



l-i 



!)• 



(40) 
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Fig. 9. Residual entropy as a function of a. The inset show a 
double logarithmic plot. 



for ai < a < aj_ x . Using cq. (40), M(C^ a) ) = 
-M(C*^' T) ), and M{C { Jj l) ) = M(C^ T) ) = 0, we obtain 

i 2 



(A|(n)} 8 



1 



2N n 



2 n-l+l 



(4 



i-i 



1) 



(43) 



In the thermodynamic limit, the correlation function of 
the Ith sector is obtained as follows: 



(A|) g = hm (A|(n)) g 



4 Z (1 



12 



(44) 



In Fig. 11, we compare the ground state correlation 
function given by eq. (44) with the low-temperature limit 
obtained by eq. (25). The agreement between both re- 
sults is reasonable. The short-range order grows expo- 
nentially as a function of sector number I. The most char- 
acteristic feature of the present system is the highly de- 
veloped short-range ordering in the intermediate phases, 
although there are no thermal and quantum fluctuations 
in the classical Ising system at T = 0. 



The magnetization per spin for the Zth sector of the nth 
stage, m n i, is written as 



M(C& T) ) 



•in— It aI—X 



(4' 



1) 



1T1-1 



(41) 



We show the a dependence of magnetizations for n = 
6, 8, 10 in Fig. 10. We have m = 1 for a < 1 and m = 
for a > 2. For the region 1 < a < 2, the magnetization is 
finite in finite size systems, but it vanishes rapidly with 
an increase in the system size. Thus, we have m = for 
1 < a < 2 in the thermodynamic limit. 
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Fig. 10. Magnetization at T = as a function of a for finite size 
systems with n = 6, 8, 10. 



5-4 Ferromagnetic correlation function 

We consider the ferromagnetic correlation function in 
the ground state of the n-stage lattice. Denoting the av- 
erage in the ground-state manifold as (■ ■ • ) g , we have 

(A|(n)) g 



^ £ [a A + a B +M{C^r B) ) 



(42) 



Fig. 11. Log-log plot of the ferromagnetic correlation function at 
T = as a function of a. 



6. Temperature Evolution of Dominant Spin 
Configurations 

Here, we study the temperature evolution of domi- 
nant spin configurations at low temperatures for the most 
characteristic case of a — 1 by considering the tempera- 
ture dependence of the entropy and the correlation func- 
tion. 

In the previous section, we have discussed the origin of 
the step structure of the residual entropy for 1 < a < 2, 
and spin configurations on each of the plateaus have 
been specified. We denote the set of the lowest-energy 
spin configurations of the Ith sector in the thermody- 
namic limit, n — * 00, as C\. On the one hand, in §4, 
we have found that the low-temperature peaks in the 
specific heat for a = 1 are related to the steps of the 
residual entropy as a function of a. Therefore, we expect 
that dominant spin configurations between two succes- 
sive low-temperature peaks can be specified in terms of 
{Ci}- ^ 

In Fig. 12, we show the temperature dependence of the 
entropy S and correlation function (A|) for a = 1. We 
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Fig. 12. Temperature dependence of (a) entropy and (b) correla- 
tion function for a = 1. 

find that a step-like structure appears in both curves at 
T/J < ^0.14. In order to specify the spin configurations 
on each of the plateaus, we present contour plots of the 
entropy and correlation function in Fig. 13. As seen from 
these contour plots, the dominant spin configurations for 
0.10 < T/J < 0.14 arc those in C 5 . With a farther de- 
crease in temperature, the dominant spin configurations 
vary as C5 — > C@ — > CV — > C$ — > ■ ■ • . 

7. Conclusions 

We have proposed a frustrated Ising hierarchical model 
with a diamond-shaped basic unit comprising four ferro- 
magnetic bonds and one non-iterated diagonal bond, and 
we have obtained a rigorous partition function. From this 
partition function, we have found that this model has 
three types of ground states depending on the magni- 
tude of the frustration parameter a. The ground state is 
ferromagnetic for a < 1 and paramagnetic with a van- 
ishing ferromagnetic correlation function for a > 2. An 
intermediate phase appears for 1 < a < 2 in our ground 
state phase diagram, and the ground-state and thermal 
properties exhibit many interesting phenomena in this 
region. When a decreases from 2 to 1, the system goes 
through an infinite number of first-order transitions to 
develop a ferromagnetic short-range order: the ferromag- 
netic correlation function of the ground-state manifold 
exhibits step-like increments and diverges, and the zero- 
temperature entropy exhibits step-like decrements and 
vanishes. Because the short-range order is highly devel- 
oped in each set of the lowest-energy spin configurations 
for 1 < a < 2, we have identified the intermediate phase 



(a) „ entropy 




10" 3 Iff 2 10"' 10° 

T/J 



(b) n correlation function 
in" 




Fig. 13. Contour plots of (a) log 5 and (b) log(A|) in the a-T 
plane, where the number corresponding to each region represents 
the sector number I. 

to be classical spin- liquid. For a — > 1, there are temper- 
ature regions in which one of these sets of spin configu- 
rations becomes the dominant one. 

In conclusion, the frustrated Ising model on the di- 
amond hierarchical lattice provides a new prototypical 
example to describe frustrated phenomena. 

Appendix A: Lowest-Energy Spin Configura- 
tions and Degeneration Factor 

Here, we study Gf?"' aB ' for a general value of n (> 
3). It is convenient to write an internal configuration as 
s = (sai,sbi,cti,sa2, sb2, 02), where s v j denotes a spin 
configuration of internal spins between the contact sites 
rj (= A or B) and j (= 1 or 2). 

For the ferromagnetic sector, setting a a = &b = we 
simply obtain 

C nj = \ s\sai,sbx G C^'pax = a; 

sa2,sb2 e C^ ) f ;a 2 = a\ . (A-l) 

Next, we consider the sector I = n — 1. In this case, 
C^n-i" 1 are constructed by the (n ~ l)-stage ground 
states of the ferromagnetic and (n — 2)-th sectors. For 
&A = &b = we obtain 

C n,n-1 = \S\SA1,SB1 € C^'j, (J X = <T\ 

SA2, S B 2 S Ln-l,n-2i G 2 = -o\ 
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+ {s\s A i, Sbi G C^1^J_ 2 ; en = -cr; 

SA2,s B2 eC ( n a f 1 > J ; f j 2 = <jy (A-2) 
For a a = —<Jb = c, we have 

f I ^,(cr,cr) — cr,cr) 

\s\sai G l7„-i,/! sbi g Q_i,„_ 2 ; °i = °"; 
SA2 G Q_ liI /„ 2 ;s_B2 G C;_ li/ ; ;(T 2 = -crj 
+ |s|sai G C^^_ 2 ; sbi G C^'J^; o-i = -cr; 

SA2 G Cttf, sb2 G C^;l 2 ; *2 = «t} . (A3) 

The first and second terms in the right-hand sides of 
eqs. (A-2) and (A-3) are graphically shown in Fig. A-l. 
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Fig. A-l. Schematic representation of spin configurations in (a) 

C nn-1 and ( b ) C< nn~l- ThG 

open (closed) circles indicate that 
these spins take the value of a (— cr). The black diamond repre- 
sents an element in C^_^f or C^— l ' f^^ ' and ^ e § ra y diamond 



represents one that in C^Z\ ^—2 or ^n-l'I-S' 



We denote the number of elements in a set A as \X\ 
then, eqs. (A-2) and (A-3) result in 



\C 



\ n (a,-<?)\ _ |,-,(<7,-cr) ,2 
l°n,n-l I ~ Z l°ri-l,n-2l > 

■f(cr,cr) 



(A-4) 



where we have used the relations, |C„_1 1 j\ = 1 and 

\ciTCLl = \ctZl2V From eq. (A-4), | is 

independent of the surface spin state. Therefore, using 

\C { °*f B) \ = 2 a ", we rewrite eq. (A-4) as 



a n = 2a„_i + 1, 



which gives 



2 n-l _ 1 



(A-5) 



(A-6) 



where we have used a% = 1. 

Finally, for a sector with I = 1, • ■ ■ , n — 2, we obtain 



{ 



S | S ai g c^vr w g ct^ru =u; 



The total number of elements in this set is 

V |/'"r( 'B,T) I — 0\B ,T) I 

X l°n-1,; ll^n-l.Z I' 



(A- 



which shows that \C 



(cAjCb) I 

.2 



does not depend on 



(<7^,(7 B ). Using |^7' CTb) | = 2 b -<, eq. (A-8) gives 

b n ,i = 46„_i,i + 1- (A-9) 
We use this relation with frj+i,; = a/+i = 2' — 1 to obtain 



\n—l — l I r*l 



1 

37 "3' 

Summarizing the above results, we find that 



(A-10) 



log|C^^| 



AU—l — l 



1 



log 2 (A-ll) 



for I — 1, 2, • ■ • 



Appendix B: Residual Entropy at Phase 
Boundary 

In order to calculate S(ai), we write a set of internal 
configurations of lowest-energy states with the surface 
spin (<ja,&b) on the n-stage lattice with the diagonal 
interaction a; as Cr? A,<JB \ai). 

Assuming a = a;, we focus on the (/ + 2)-stage lat- 
tice. Because £'(+2(0;;) = E l+2 .i = E l+2 ^ +1 , we obtain 
Cll A 2 f B) C Cfe'™>(aO and c£## C C^VO- 
Note that the set C^^'i is constructed from four ele- 



ments in X)<r o-' ^1+11 1 ano ^ another set C^Z^'m! IS con ~ 
structed from two elements in CT , c[^ii and two ele- 
ments in £ cfc" } 



Now, we consider AC^ B \ai) = C^' aB> '{on) - 

^+2,f B) - C&i+i 1 ' Because £,, +a ,i = £z+2,;+i at a = 
011, we notice that there are other lowest-energy states 
with three elements in ^2 aa i ^j+iV and one clement in 

-i(cr,cr') 



AC, 



{<J A ,U B ) 

1+2 



("0 



Jolt, r- /~<(<* A ,0 a) . „ f- s~l(,0B,crA) . „ „ . 

\s\sa\ G C 1+1 j ,s B1 €C |+lil ;ax = CTA\ 

SA2 G C l+l l ; s B2 G C l+l l ; cr 2 = -cr^j 

ill r- /^(^A^b) ~- s~l(o B ,cr b) _ _ 

+ |s|sai G C, v +li , ,SbiGC; v +1 j, y ;cri=cr B ; 

sa2 g ;s B2 eo ;+lj , ;cr 2 = -cr B | 

+ | s|sai G , sbi G , <xi = -cta, 

SA2 G i SB2 G + 1 >; , (7 2 = CTA J 

ill r^(. a a-~&b) ~- ri( G B-—GB) _ _ 

+ |s|sai G C; +1 j ;ssiG6f +lji y ;cri = -o- B ; 

SA2 G C, v +1:i ;SB2GC, v +1j ,02 = t7 B j .(B-l) 

Schematic representations of the elements in the first, 
second, third, and fourth sets in the right-hand side of 
eq. (B-l) are respectively shown in Figs. B-l(a)-(d). 
Equation (B-l) gives 



AClir B \ ai )\ = 2 



3-2'-l 



(B-2) 
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Thus, the entropy at a = on is given by 
S(ai) = lim -Llog4-2 d " 

n— >oo l\ n 

= Sooi + ^ log [l + 4(2~ 2 ' + 4~ 2 ') 



Fig. B-l. Schematic representation of spin configurations in 
ACj^2 ' <7B \cti). Open and closed circles indicate that the di- 
rections of those spins are opposite to each other. 



and thus, we obtain 

\C\l A 2 aB \ ai )\ = 2 4 ' 2 '- 3 + 2 2 ' +1 - 1 + 2 3 ' 2 '- 1 . (B-3) 
For the n-stage lattice with n > I + 2, using 

\C&w){a l )\=2 d ~, (B-4) 

we obtain 

d n = 4d„_! + 1. (B-5) 
The recursion relation with the initial condition 

^ = 4.2'-3 + '^ 1 + 2 rr +2 "' > (B-6) 
gives 



log 2 



, . 8 log(l + 2 
4-2' i — 



-2 l + 1 _|_ 2 2 -2 ! 



log 2 



1 

(B-7) 
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